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Abstract. Quantum computation and quantum information are of great 
current interest in computer science, mathematics, physical sciences 
and engineering. They will likely lead to a new wave of technological 
innovations in communication, computation and cryptography. As the 
theory of quantum physics is fundamentally stochastic, randomness and 
uncertainty are deeply rooted in quantum computation, quantum simu- 
lation and quantum information. Consequently quantum algorithms are 
random in nature, and quantum simulation utilizes Monte Carlo tech- 
niques extensively. Thus statistics can play an important role in quan- 
tum computation and quantum simulation, which in turn offer great 
potential to revolutionize computational statistics. While only pseudo- 
random numbers can be generated by classical computers, quantum 
computers are able to produce genuine random numbers; quantum 
computers can exponentially or quadratically speed up median eval- 
uation, Monte Carlo integration and Markov chain simulation. This 
paper gives a brief review on quantum computation, quantum simu- 
lation and quantum information. We introduce the basic concepts of 
quantum computation and quantum simulation and present quantum 
algorithms that are known to be much faster than the available clas- 
sic algorithms. We provide a statistical framework for the analysis of 
quantum algorithms and quantum simulation. 

Key words and phrases: Quantum algorithm, quantum bit (qubit), 
quantum Fourier transform, quantum information, quantum mechan- 
ics, quantum Monte Carlo, quantum probability, quantum simulation, 
quantum statistics. 



1. INTRODUCTION 

For decades computer hardware has grown in pow- 
er approximately according to Moore's law, which 
states that the computer power doubles for con- 
stant cost roughly once every two years. However, 
because of the fundamental difficulties of size in con- 
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ventional computer technology, this dream run is 
ending. The conventional approaches to the fabrica- 
tion of computer technology are to make electronic 
devices smaller and smaller in order to increase the 
computer power. As the sizes of the electronic de- 
vices get close to the atomic scale, quantum effects 
are starting to interfere in their functioning, and 
thus the conventional approaches run up against the 
size limit. One possible way to get around the dif- 
ficulties is to move to a new computing paradigm 
provided by quantum information science. Quantum 
information science is based on the idea of using 
quantum devices to perform computation and ma- 
nipulate and transmit information, instead of elec- 
tronic devices following the laws of classical physics, 
see Deutsch (1985), DiVincenzo (1995), Feynman 
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(1981/82). Quantum mechanics and information the- 
ory are two of the great scientific developments and 
technological revolutions in the 20th century, and 
quantum information science is to marry the two 
previously disparate fields and form a single unify- 
ing viewpoint. Quantum information science studies 
the preparation and control of the quantum states 
of physical systems for the purposes of information 
transmission and manipulation. It includes quantum 
computation, quantum communication and quan- 
tum cryptography. This revolutionary field will en- 
able a range of exotic new devices to be possible. 
There is now a general agreement that quantum in- 
formation science will likely lead to the creation of 
a quantum computer to solve problems that could 
not be efficiently solved on a classical computer. 
Already scientists have built rudimentary quantum 
computers in the research laboratory to run quan- 
tum algorithms and perform certain calculations. In- 
tensive research efforts are under way around the 
world to investigate a number of technologies that 
could lead to more powerful and more prevalent 
quantum computers in the near future. It is believed 
that quantum information and quantum bits are to 
lead to a 21st century technological revolution much 
as classic information and classic bits did to the 
20th century. Since the theory of quantum mechan- 
ics is fundamentally stochastic, randomness and un- 
certainty are deeply rooted in quantum computation 
and quantum information. As a result, quantum al- 
gorithms are of random nature in the sense that 
they yield correct solutions only with some prob- 
abilities, and Monte Carlo methods are widely em- 
ployed in quantum simulation. Thus statistics has 
an important role to play in quantum computation, 
quantum simulation and quantum information. On 
the other hand, quantum computation and quantum 
simulation have tremendous potential to revolution- 
ize computational statistics. 

A quantum system is generally described by its 
state, and the state is mathematically defined to be 
a unit vector in some complex Hilbert space. The 
number of complex numbers required to characterize 
the quantum state usually grows exponentially with 
the size of the system, rather than linearly, as occurs 
in classical systems. As a consequence, it takes an 
exponential number of bits of memory on a classical 
computer to store the quantum state, which puts 
classical computers in a difficult position to simu- 
late a quantum system. On the other hand, nature 
quantum systems are able to store and keep track 
of an exponential number of complex numbers and 



perform data manipulations and calculations as the 
systems evolve. Quantum information science is to 
grapple with understanding how to take advantage 
of the enormous information hidden in the quan- 
tum systems and to harness the immense potential 
computational power of atoms and molecules for the 
purpose of performing computation and processing 
information. Already it has been shown that quan- 
tum algorithms like Grover's search algorithm and 
Shor's factoring algorithm provide great advantage 
over known classical algorithms. 

Contemporary scientific studies often rely on un- 
derstanding complex quantum systems, such as those 
in biochemistry and nanotechnology for the design 
of biomolecules and nano- materials. Quantum sim- 
ulation is to use computers to simulate a quantum 
system and its time evolution. Classical computers 
are being used for quantum simulation in design- 
ing novel molecules and creating innovative nano- 
products. Quantum computers built upon quantum 
systems may excel in simulating naturally occurring 
quantum systems, while large quantum systems may 
be impossible to simulate in an efficient manner by 
classical computers. A quantum system with b dis- 
tinct components may be described with b quantum 
bits in a quantum computer, while a classical com- 
puter requires 2^ bits of memory to store its quan- 
tum state. This advantage allows quantum comput- 
ers to efficiently simulate general quantum systems 
that are not efficiently simulatable on classical com- 
puters. 

In this article we review the concepts of quantum 
computation and introduce quantum algorithms and 
quantum simulation. The quantum algorithms are 
known to be much faster than the available classi- 
cal algorithms. Statistical analyses of quantum al- 
gorithms and quantum simulation are provided. We 
give a brief description on quantum information. 
The article sections start with presentations in broad 
brushstrokes, followed by specific discussions along 
with some mathematical derivations if necessary. The 
intention is to give each topic first an overview and 
then a general description and a precise characteri- 
zation. It is recommended to focus on the qualitative 
discussions but skip the derivations for the readers 
who would like to get a quick picture of quantum 
computation and quantum simulation. 

The rest of the paper proceeds as follows. Sec- 
tion 2 briefly introduces quantum mechanics, quan- 
tum probability and quantum statistics. Section 3 
reviews basic concepts of quantum computation and 
entanglement. Section 4 illustrates some widely known 
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quantum algorithms and provides a statistical frame- 
work for the study of quantum algorithms. Section 5 
presents quantum simulation and discusses its sta- 
tistical analysis. Section 6 gives a short description 
on quantum information theory. Section 7 features 
concluding remarks and lists some open research 
problems. 

2. BRIEF BACKGROUND REVIEW ON 
QUANTUM THEORY 

Quantum mechanics has been applied to every- 
thing under and inside the Sun, from chemical reac- 
tion and superconductor to the structure of DNA 
and nuclear fusion in stars. Although the signifi- 
cant difference between classical physics and quan- 
tum physics lies in the quantum prediction of phys- 
ical entity when the scale of observations becomes 
comparable to the atomic or sub-atomic scale, many 
macroscopic properties of systems can only be fully 
explained and understood by quantum physics. The 
quantum world is extremely strange, and quantum 
theory is completely counterintuitive. Light waves 
behave like particles and particles behave like waves 
(wave particle duality) ; matter can go from one spot 
to another without moving through the intermedi- 
ate space (quantum tunneling); information can be 
moved across a vast distance without transmitting 
it through the intervening space (quantum telepor- 
tation). Quantum theory provides a mathematical 
description of wave particle duality and interaction 
of matter and energy. It describes the time evolu- 
tions of physical systems via wave functions. The 
wave functions encapsulate the probabilities that 
particles are to be found in a given state at a given 
time. For example, the probability of finding a pho- 
ton in some region is the square of the modulus 
of a wave function, and, since at some point the 
sum of two wave functions can be zero but neither 
wave function is zero, probabilities appear to can- 
cel out each other in a way totally unexpected from 
classical probability. The intrinsic stochastic nature 
of quantum theory indicates a deep connection be- 
tween quantum mechanics and probability. Since the 
main focus of this paper is on quantum computation 
and quantum information, we give a brief descrip- 
tion of quantum theory in this section to provide 
some quantum background for the purpose of re- 
viewing quantum computation and quantum sim- 
ulation in subsequent sections. For further reading 
on the subjects we recommend textbooks by Saku- 
rai and Napolitano (2010) at the graduate level and 
Griffiths (2004) at the undergraduate level for quan- 



tum mechanics, Holevo (1982), Parthasarathy (1992) 
and Wang (1994) for quantum probability and quan- 
tum stochastic processes, and Artiles, Gill and Guta 
(2005) and Barndorff-Nielsen, Gill and Jupp (2003) 
for quantum statistics. 

2.1 Hilbert Space and Operator 

For the sake of simplicity we choose to work with 
comparatively easy finite-dimensional situations. De- 
note by C the set of all complex numbers. We start 
with vector space in linear algebra. A simple exam- 
ple of vector space is C'"' consisting of all fe-tuples 
of complex numbers (zi, . . . , z^). The elements of 
a vector space are called vectors. As in quantum 
mechanics and quantum computation, we use Dirac 
notations |-) (which is called ket) and (-1 (which is 
called bra) to indicate that the objects are column 
vectors or row vectors in the vector space, respec- 
tively. Denote by superscripts *, / and f the conju- 
gate of a complex number, the transpose of a vec- 
tor or matrix, and conjugate transpose operation, 
respectively. We define an inner product on the vec- 
tor space to be a function that takes as input two 
vectors from the vector space and produces a com- 
plex number as output. For |u) and \v) in the vec- 
tor space, we denote their inner product by {u\v). 
The inner product must satisfy (i) conjugate sym- 
metry, {u\v) = {{v\u))*; (ii) linearity in the second 
argument, {u\v + tt;) = {u\v) + {u\w); (iii) positive- 
definiteness, {u\u) > with equality only for u = 0. 
For example, C'' has a natural inner product 

k 

{u\v) = ^ U*Vj = {ul,...,ul){vi,..., Vk)\ 

where (m| = (ui, . . . , Ufc) and = (ui, . . . , Vf.)' . An 
inner product induces a norm \\u\\ = \J and 
a distance ||n — v\\ between |n) and \v). For the 
finite-dimensional case, a Hilbert space is simply 
a vector space with an inner product. 

An operator A on 7^, denoted by A(|ii)) for \u) G 7^, 
is a function mapping from Ti to Ti that satisfies 
A{a\u) + b\v)) = aA(|'u)) + 6A(|t>)) for any \u), \v) G V. 
and a,b gC We can represent an operator through 
a matrix. Suppose that A is an operator on Ti and 
ei , . . . , efc form an orthonormal basis in H . Then 
each A(|ej)) G Ti and there exists a unique k x k 

matrix [ajc) such that A(|ej)) = X^^^^ \ei)aij, j = 
l,...,k. We will identify operator A with matrix 
(aji) and use A for both operator and matrix (aji). 
An operator A on ^ is said to be self-adjoint if 
its corresponding matrix A is Hermitian, that is. 
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A = A^^. We also refer to self-adjoint operators as 
Hermitian operators. An operator U is said to be 
unitary if its corresponding matrix U is unitary, 
that is, UlJl^ = TJl^U = I. We say an operator A is 
semi-positive (or positive) definite if its correspond- 
ing matrix A is semi-positive (or positive) definite, 
that is, {u\p\u) > for \u) G H (or {u\p\u) > for 
\u) ^Ti with equality only for |n) =0). The trace 
of an operator A, denoted by Tr(A), is defined to 
be the trace of its corresponding matrix A = {aji), 
that is, Tr(A) = Yl^=i '^jj- 
2.2 Quantum System 

Quantum mechanics depicts phenomena at micro- 
scopic level such as position and momentum of an in- 
dividual particle like an atom or electron, spin of an 
electron, detection of light photons, and the emission 
and absorption of light by atoms. Unlike classical 
mechanics where physical entities like position and 
momentum can be measured precisely, the theory of 
quantum mechanics is intrinsically stochastic in a 
sense that we can only make probabilistic prediction 
about the results of the measurements performed. 

Quantum mechanics is mathematically described 
by a Hilbert space H and self-adjoint operators onH. 
A quantum system is completely characterized by its 
state and the time evolution of the state. A state is 
defined to be a unit vector in %. Let \ip{t)) be the 
state of the quantum system at time t, which is also 
referred to as a wave function. The states \tp{ti)) and 
\ip{t2)) at ti and t2 are connected through |V'(i2)) = 
U(ti, t2)|V'(^i))5 where U(ti,t2) is a unitary opera- 
tor depending only on time ti and t2- In fact, there 
exists a self-adjoint operator H, which is known as 
the Hamiltonian of the quantum system, such that 
U(ti,t2) = exp[— iH(t2 — ti)]- With Hamiltonian H, 
we may describe the continuous time evolution of 
\ilj{t)) by the Schrodinger equation 

(1) i3|M = HK.W). i = v/=T, 

Alternatively a quantum system can be described 
by a density operator (or density matrix). A density 
operator p is an operator on % which (1) is self- 
adjoint; (2) is semi-positive definite; (3) has unit 
trace [i.e., Tr(p) = 1]. Following the convention in 
quantum information science, we reserve notation p 
for state, density operator or density matrix. A state 
is often classified as a pure state or an ensemble of 
pure states. A pure state is a unit vector jV') in ^, 
which corresponds to a density operator p= 
and an ensemble of pure states corresponds to the 



case that the quantum system is in one of states 
iV'j), J = 1, . . . , J, with probability pj being in state 
I'f/'j), and the corresponding density operator 

J 

(2) p=Y,vMj){n- 

See Griffiths (2004), Sakurai and Napolitano (2010) 
and Shankar (1994). 

2.3 Quantum Probability 

We can test the theory of quantum mechanics by 
checking its predictions with experiments of per- 
forming measurements on quantum systems in the 
laboratory. The usual quantum measurements are 
on observables such as position, momentum, spin, 
and so on, where an observable X is defined as a self- 
adjoint operator on Hilbert space %. The observable 
definition is motivated from the fact that the eigen- 
values of self-adjoint operators are real. Assume that 
an observable X has a discrete spectrum with the 
following diagonal form 

p 

(3) X = ^XaQa, 

a=l 

where Xa are real eigenvalues of X and Qa are the 
corresponding one-dimensional projections onto the 
orthogonal eigenvectors of X. Consider such an ob- 
servable in the quantum system prepared in state p. 
Measure space iyt,F) is used to describe possible 
measurement outcomes of the observable, and the 
result of the measurement is a random variable on 
(O, J^) with probability distribution Pp. We denote 
by X the result of the measurement of observable X 
given by (3). Then X is a random variable taking 
values in {xi, a;2, ...,}, and under pure state IV'), the 
probability that measurement outcome Xa occurs is 
defined to be 

P{a) = Pp{X = Xa) 

= (V'IQa|V')=Tr(Q,|^)(^|), a = 1,2,.... 

With the probability we derive the expectation un- 
der pure state IV') , 

E^(X) = Y,XaP{a) = Y,Xa{^\W) 

a a 

= (V|X|V^)=TV(X|V)(Vl). 

Note the difference between an observable X which 
is a Hermitian matrix and its measurement result X 
which is a real-valued random variable. 

Measuring observable X will alter the state of 
the quantum system (Kiefer (2004); von Neumann 
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(1955)). If the quantum system is prepared with ini- 
tial state I , the state of the system after the mea- 
surement result Xa is defined to be 



(4) 



Qa\ 



For an ensemble state with density operator p 
given by (2), if the quantum state is {tpj), the prob- 
ability that result Xa occurs is 

P{a\j) = (VilQalV',) = Tr(Q„|V',)(V^,|). 

Applying the law of total probability, we obtain that 
under state p, the probability that Xa occurs is equal 
to 



P{a) = Pp{X = Xa) 



J 

E 



PjP{a\j) 



= ^p,TY(QJVi>(V'il) = Tr(Q,p). 
j=i 

The expectation of X under state p, 

p p 
Ep[X] = XaPp[X = Xa]=Y, Xa Tr(Qap) 

a=l a=l 

= tr(Xp), 

and variance 

Varp[X]=tr[X2p]-(tr[Xp])2. 

We may derive the density operator of the quan- 
tum system after obtaining the measurement re- 
sult Xa by conditional probability arguments as fol- 
lows. If the quantum system is in pure state 
before the measurement, the quantum state after 
measurement result Xa has occurred is 

If the quantum state is p before the measurement, 
after observing measurement outcome Xa we have 
the following ensemble of states: the quantum sys- 
tem is in pure state IV'") with probability P(j|a), 
where Bayes's theorem shows 

P{j\a)=pjP{a\j)/P{a). 

Thus after measurement Xa the density operator for 
the ensemble state is given by 

J 

Pa = Y.p{3\a)\rj){r^ 



Qa|V'j)(V'ilQfi 



Qa\lpj){'>Pj\Qa _ QaPQa 



P{a) 



Tr(Q„p)- 



See Holevo (1982), Parthasarathy (1992) and Saku- 
rai and Napolitano (2010). 

2.4 Quantum Statistics 

For a given quantum system, it is very important 
but difficult to know its state. If we do not know in 
advance the state of the quantum system, we may 
infer the quantum state by the measurement results 
of some observables obtained from the quantum sys- 
tem and show that a certain state has been created. 
In statistical terminology, we want to estimate den- 
sity matrix p based on measurements on an often 
large number of systems which are identically pre- 
pared in the state p. That is, after measuring observ- 
ables on some identical quantum systems, we can 
make statistical inference about probability distri- 
bution Pp of the measurements and thus indirectly 
about density matrix p. In the literature of quantum 
physics, quantum tomography is referred to as the 
reconstruction of the underlying density matrix p by 
probing identically prepared quantum systems from 
some different angles. Specifically, suppose that we 
perform measurements of observables on identically 
prepared quantum systems in an unknown state p 
and obtain measurement results Xi , X^- Assume 
that p is known up to some unknown parameter 9; 
then Xi, . . . ,Xn are i.i.d. observations with distri- 
butions Pp which depend on 9. This gives a quan- 
tum parametric statistical model. We may then de- 
fine quantum likelihood and Fisher quantum infor- 
mation and establish quantum point estimation and 
quantum hypothesis testing theory. Alternatively we 
may model p nonparametrically by assuming that p 
is an infinite matrix and then use nonparametric 
methods to estimate the density matrix. For details 
see Artiles, Gill and Gu^a (2005), Barndorff-Nielsen, 
Gill and Jupp (2003), Butucea, Gu^a and Artiles 
(2007) and Nussbaum and Szkola (2009). 

3. QUANTUM COMPUTING CONCEPTS 

Unlike classical computers using transistors to 
crunch the ones and zeroes individually, quantum 
computers can handle both one and zero simulta- 
neously via what are known as superposition quan- 
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turn states. A superposition state is a state of mat- 
ter which we may think of as both one and zero at 
the same time. Quantum computers use the strange 
superposition states and quantum entanglements to 
do the trick of performing simultaneous calculations 
and extracting the calculated results. The spooky 
phenomena of quantum entanglement and superpo- 
sition are the key that enables quantum computers 
to be superfast and vastly outperform classical com- 
puters. 

3.1 Quantum Bit 

Analogous to the fundamental concept of bit in 
classical computation and classical information, we 
have its counterpart, quantum bit, in quantum com- 
putation and quantum information. Quantum bit 
is called qubit for short. Just like a classical bit 
with state either or 1, a qubit has states |0) and 
However, the real difference between a bit and 
a qubit is that besides states |0) and |1), a qubit 
may take the superposition states, 

l-f/-) =ao|0) +ai|l), 

where ao and ai are complex numbers and called 
amplitudes satisfying |aoP + l^ip = 1. That is, the 
states of a qubit are unit vectors in a two-dimensional 
complex vector space, and states |0) and |1) consist 
of an orthonormal basis for the space and are often 
referred to as computational basis states. For a clas- 
sical bit we can examine it to determine whether it 
is in the state or 1. However, for a qubit we cannot 
determine its state and find the values of ao and ai 
by examining it. The stochastic nature of quantum 
theory shows that we can measure a qubit and ob- 
tain either the result 0, with probability |aoP, or 
the result 1, with probability |aip. Physical exper- 
iments have realized qubits as physical objects in 
different physical systems, such as the two states 
of an electron orbiting a single atom, the two dif- 
ferent polarizations of a photon, or the alignment 
of a nuclear spin in a uniform magnetic field. Con- 
sider the case of atom model by corresponding |0) 
and |1) with the so-called "ground" and "excited" 
states of the electron, respectively. As the atom is 
shined by light with suitable energy and for a proper 
amount of time, the electron can be moved from the 
|0) state to the |1) state and vice versa. Further- 
more, by shortening the length of time shining the 
light on the atom, we may move the electron ini- 
tially in the state |0) to "halfway" between |0) and 
say, into a state (|0) + |l))/\/2. 
Note that for qubit state the only measurable 
quantities are the probabilities |aoP and |aip; since 



|e*^aa;p = |aa:P) where x = 0,1, i = and 6 is 

a real number, from the viewpoint of the qubit mea- 
surements, states e^^l'i/') and IV') are identical. That 
is, multiplying a qubit state by a global phase fac- 
tor e*^ bears no observational consequence. 

Note the distinction between superposition states 
and probability mixtures (or ensemble of pure states 
defined in Section 2.1). Consider superposition (|0) -|- 
|l))/\/2 as a pure state. Its density matrix is given 
by 

i(|0) + |l))((0| + (l|) 

= i(|0)(0| + |l)(l|) + i(|0)(l| + |l)(0|), 

while the first term on the right-hand side of the 
above equation corresponds to the ensemble of pure 
states |0) and that is, a probabilistic mixture of 
states |0) and |1) with equal probability. 

Similar to classic bits, we can define multiple qu- 
bits. The states of b qubits are unit vectors in a 2^- 
dimensional complex vector space with 2^* computa- 
tional basis states of the form \x1X2 ■ ■ -Xb), Xj = 
or 1, j = 1, . . . ,b. For example, the states of two 
qubits are unit vectors in a four-dimensional com- 
plex vector space, with four computational basis 
states labeled by |00), |01), |10) and The com- 
putational basis states |00), |01), |10) and |11) gen- 
erate the four-dimensional complex vector space, and 
the superposition states are all unit vector in the 
space with the forms 

1^) = aoolOO) + aoilOl) + aio|10) + aii|ll), 

where amplitudes are complex numbers satisfy- 
ing |aooP + |aoiP + |aioP + |ttiiP = 1- As in the 
single qubit case, when two qubits are measured we 
get result x being one of 00,01,10,11, with prob- 
ability |q!xP- Moreover, we may measure just the 
first qubit of the two-qubit system and obtain ei- 
ther the result 0, with probability |aooP + |"oiP) 
or the result 1, with probability |aioP + As 
quantum measuring changes the quantum state, if 
the measurement result on the first qubit is 0, after 
the measurement the qubits are in the state 

aoolOO) +aoi|01) 
/I 12^1 

A qubit is the simplest quantum system. The quan- 
tum system of b qubits is described by a 2''-dimension- 
al complex vector space with each superposition state 
specified by 2^ amplitudes. As 2'' increases exponen- 
tially in b, it is very easy for such a system to have 
an enormously big vector space. A quantum sys- 
tem with even a few dozens of "qubits" will strain 
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the resources of even the largest supercomputers. 
Consider a quantum system of 50 qubits. 2^^^ « 10^^ 
complex amplitudes are needed to specify its quan- 
tum states. With 128 bits of precision, it requires 
approximately 32 thousand terabytes of information 
to store all 10^^ complex amplitudes. Such storage 
capacity may be available in future supercomputers. 
For a quantum system with b = 500 qubits we need 
to specify 2^"'' complex amplitudes for its states. It 
is unimaginable to store all 2^^*^ complex numbers 
in any classical computers. In principle, a quantum 
system with only a few hundred atoms can manage 
such an enormous amount of data and execute cal- 
culations as the system evolves. Quantum computa- 
tion and quantum information are to find ways to 
utilize the immense potential computational power 
in quantum systems. 

3.2 Quantum Circuit Model 

As a classical computer is built from an electri- 
cal circuit consisting of wires for carrying informa- 
tion around the circuit and logic gates for perform- 
ing simple computational tasks, a quantum com- 
puter can be created from a quantum circuit with 
quantum gates to perform quantum computation 
and manipulate quantum information. A number of 
physical systems are being investigated for build- 
ing quantum computers. These include optical pho- 
ton, optical cavity quantum electrodynamics, ion 
traps, nuclear magnetic resonance with molecules, 
quantum dots, and superconductors (Nielsen and 
Chuang (2000)). In fact, primitive solid-state quan- 
tum processors have been created in research lab- 
oratories to run quantum algorithms (DiCarlo et 
al. (2009); Johnson et al. (2011); Mariantoni et al. 
(2011); Sayrin et al. (2011)). The circuit model is 
particularly important in quantum computation and 
quantum information, and a quantum computer is 
often synonymous with the quantum circuit model. 
A quantum circuit operates on b qubits for some in- 
teger b. The state takes a form of \xi---Xb), with 
state space being a 2^-dimensional complex Hilbert 
space. When Xj = or 1, states \xi---Xb) are the 
computational basis states of the quantum computer 
and often written as \x), where x is the integer with 
binary representation xi - ■ -Xb- 

As a classical logic gate converts classical bits from 
one form to another such as — )• 1 and 1 — )• 0, a quan- 
tum gate operates on qubits. Quantum mechanics 
dictates that quantum gates operating on b qubits 
are 2^ by 2^ unitary matrices on the 2''-dimensional 
Hilbert space. For example, a Hadamard gate on one 



qubit is the 2x2 unitary matrix that realizes the 
following transformation: 



|0) 



|0) + |1) 



ID 



in 



V2 ' ' ' ^/2 ■ 
Consider another important gate on two qubits which 
is called control-NOT gate. It takes the two input 
qubits as control qubit and target qubit, respec- 
tively, and the output target qubit of the gate re- 
tains the input target qubit if the control qubit is 
|0) and is flipped if the control qubit is |1), that is, 

|00)^|00), |01)^|01), 

|10)^|11), |11)^|10). 

Generally for any single qubit unitary operation U, 
a control-U gate is a two-qubit gate, with one con- 
trol qubit and one target qubit. If the control qubit 
is |1), U is applied to the target qubit; if the control 
qubit is |0), the target qubit is left alone, that is, 

|0)|0)^|0)|0), |0)|1)^|0)|1), 
|1)|0)^|1)U|0), |1)|1)^|1)U|1). 

If f{x) maps {0, 1}'' onto {0, 1}, we define a unitary 
transformation Uj that operates on 1 qubit state 

(6) \x,y) ^\x,y® f{x)), 

where x = xi - ■ ■ Xf, with xj = or 1 is the data reg- 
ister, y = or 1 is the target register, © denotes ad- 
ditional modulo 2. If ?/ = 0, after the transformation 
Uj, the state of the last qubit is the value of f{x). 

3.3 Entanglement 

Quantum entanglement is one of the most mind- 
bending creatures known to science. It is referred 
to as the phenomenon that two qubits behave like 
twins that are connected by an invisible wave to 
share each other's properties. 

3.3.1 Bell states Consider a quantum gate on two- 
qubit basis states |00), |01), |10) and |11) that is 
composed of a Hadamard gate on the first qubit and 
then is followed by a control-NOT gate. The output 
states of the gate are as follows: 



|00) 



1 10) 



|00) + |11) 

V2 
|00)-|11) 



|01) 



11 



|01) + |10) 

V2 
|01)-|10) 



V2 ' ' ' V2 ' 

Physicists Bell, Einstein, Podolsky and Rosen dis- 
covered the amazing properties of these four states, 
which are often referred to as the Bell states, EPR 
states or EPR pairs (Bell (1964); Einstein, Podolsky 
and Rosen (1935)). In general states such as these 
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four states that cannot be expressed as products of 
some single qubits are called entangled states. En- 
tangled states, which are not fully understood in 
quantum physics, have remarkable properties. 
For the two-qubit system consider a Bell state 

i.a |oi)-|io) 



and an observable 



M = UxO-^ + ay(Ty + az(Tz, 

where {ax,ay, a^) is a real unit vector (i.e., + + 
= 1), and <t^., cTy and are Pauli matrices given 

by 



(7) 



-1 



-1. 



It is easy to show that M has eigenvalues ±1 for any 
real unit vector {ax,ay,az)- If measuring observable 
M on each qubit of {tp), we will obtain a measure- 
ment result of +1 or —1. Surprisingly, no matter 
what choice of (ox, Oy, o^), the measurement results 
on the two qubits are always opposite of each other, 
that is, when the first qubit measurement is —1, then 
the second qubit measurement will be +1, and vice 
versa. 

The two-qubit system can be realized by the spins 
of two particles, and the measurement of M is re- 
ferred to as a measurement of spin along the {ax, ay, 
az) axis. After the two-particle system is prepared in 
the Bell state {ip), the two particles drift far apart. 
Alice measures the spin of the first particle and Bob 
measures the spin of the second particle. The above 
opposite measurement phenomenon corresponds to 
that due to the entangled state if Alice gets 
a result +1 from her spin measurement on the first 
particle, then the state of the system immediately 
jumps to the untangled state so that the second par- 
ticle now has definite spin state and Bob's spin mea- 
surement on the second particle gives definite result 
— 1. This phenomenon is often referred to as anti- 
correlation in entanglement experiments (Neumann 
et al. (2008); Sakurai and Napolitano (2010)). 

The mathematical arguments for the anti-correla- 
tion phenomenon are as follows. The measurement 
of M on the first (or second) qubit of lip) corre- 
sponds to the spin measurement of Alice's (or Bob's) 
particle along the {ax,ay,az) axis in the above two- 
particle spin model. From Sections 2.3 and 3.1 we 



have that the two-qubit system is described by the 
Bell state in C^; measuring M on the first qubit 
of means performing measurement on observable 
M (8) I in the Bell state \ip), which alters the quan- 
tum state of the two-qubit system; measuring M on 
the second qubit corresponds to measuring observ- 
able I (8) M in the altered quantum state, where I is 
the 2 by 2 identity matrix, and M (g) I and I (8 M 
are matrix tensor products. 

Denote by |<^±) the two orthonormal eigenvectors 
of M corresponding to eigenvalues ±1, respectively, 
and let Q-t be the respective projections onto the 
eigenvectors \(p±). Following (3)-(5) we have a diag- 
onal representation M = Q+ — Q_; when we mea- 
sure observable M on each qubit, the possible mea- 
surement results are ±1; measuring M on the first 
qubit changes the state of the two-qubit system, and 
after the measurement result ±1 on the first qubit, 
the post-measurement state of the two-qubit sys- 
tem is Q±'8'I|'!/')/IIQ± Below we will evalu- 
ate the post-measurement state and show that mea- 
suring I (8) M in the post-measurement state always 
yields measurement results opposite to the measure- 
ment results on the first qubit. 

Since (|0),|1)) and (I93+), |(/5_)) are two bases for 
the one-qubit system in C^, then 



|0) 

ID 



an ai2 
021 022 



where {oiji) forms a 2 x 2 unitary matrix with de- 
terminant equal to a phase factor e*^ {i = \J—1) for 
some real 6. Substituting the above expressions into 
the entangled state and ignoring a global phase fac- 
tor e*^ (which has no effects on measurement results; 
see Section 3.1), we obtain 

|01)-|10) 



(8) 



V2 



V2 



V2 



From the definitions oi\(p±) and Q-t , Q+ (8 1 1 9?- ) = 

0, and Q_ (8l|(/3+(/9_) = 0. If the measurement result 
of M on the first qubit is +1 (or —1), from (8) we 
obtain the post-measurement state of the two-qubit 
system as follows: 



Q 



Q+ (8 - Q+ (8 l\(p-ipj, 



IQ- 



V2IIQ- 
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(or Q_®I|V^)/||Q_®I|V^)|| =6^^-^^+) ~ 
Since 

= -Iv'+v^-), 

that is, the post-measurement state (or 
is the eigenvector of I (8) M corresponding 
to eigenvalue —1 (or +1), performing measurement 
on I (g) M in the post-measurement state must al- 
ways yield measurement result —1 (or -fl). Thus, 
the measurement results of M on the two qubits of 
are always opposite to each other. 

3.3.2 Quantum teleportation Quantum teleporta- 
tion is a process by which we can transfer the state 
of a qubit from one location to another, without 
transmitting it through the intervening space. We 
illustrate the phenomenon as follows. Alice and Bob 
together generated a Bell state long ago. Each took 
one qubit of the Bell state when they split. Now they 
are far away from each other. The mission for Alice is 
to deliver a qubit |^) to Bob, while he is hiding, and 
she can only send classical information to Bob but 
does not know the state of the qubit \^). Quantum 
teleportation is a way that Alice utilizes the entan- 
gled Bell state to send a qubit of unknown state to 
Bob, with only a small overhead of classical com- 
munication. Recently a breakthrough in quantum 
teleportation has been made by successfully trans- 
ferring complex quantum data instantaneously from 
one place to another, paving the way for real- world 
applications of quantum communications (Lee et al. 
(2011)). 

Here is how it works. Alice interacts the qubit {ip) 
to be teleported with her half of the Bell state, and 
then performs a measurement on the two interacted 
qubits to obtain one of four possible two-classical- 
bit results: 00,01,10 and 11. She sends the two-bit 
information via classical communication to Bob. De- 
pending on Alice's classical message. Bob performs 
one of four operations on his half of the Bell state. 
Surprisingly, the described procedure allows Bob to 
recover the original state \ip). 

Specifically assume that the state to be teleported 
is {ip) = ao|0) -|-ai|l), where qq and ai are unknown 
amplitudes. First, consider a three-qubit state 

I , J00) + |11) 
= 1^/ ^ 

= i=[ao|0)(|00) + 111)) +ai|l)(|00) + |11))], 



where the first two qubits (on the left) belong to 
Alice, and the third qubit to Bob. Note that Alice's 
second qubit and Bob's third qubit are from the 
entangled Bell state. Second, Alice applies a control- 
NOT gate to her qubits in {ipo) and obtains 

= 7i[«o|0)(|00) + 111)) +ai|l)(|10) + |01))]. 

Third, she applies a Hadamard gate to the first qubit 
in \ipi) and gets 

|^2) = i[ao(|0) + |l))(|00) + |ll)) 

+ ai(|0)-|l))(|10) + |01))]. 

We regroup the terms of \<p2) and rewrite it as fol- 
lows: 

\^2) = i[|00)(ao|0) +ai|l)) + |01)(ao|l) + ai|0)) 

+ |10)(ao|0)-ai|l)) 

+ |ll)(ao|l)-ai|0))]. 

The new expression has four terms, and each term 
has Alice's qubits in one of four possible states |00), 
|01),|10) and |11), and Bob's qubit is in the state 
related to the original state \tl^). If Alice performs 
a measurement on her qubits and informs Bob of 
the measurement result, then his post-measurement 
state is completely determined. For example, the 
first term has Alice's qubits in the state |00) and 
Bob's qubit in state IV'). Therefore, if Alice's mea- 
surement result on her qubits is 00, then Bob's qubit 
will be in state {ip). Below is a list of Bob's four post- 
measurement states corresponding to the results of 
Alice's measurements: 

OO^aolO) -Fai|l), 01 ^ ao|l) + ai|0), 

lO^aolO) -ai|l), ll^ao|l)-ai|0). 

As Alice's measurement outcome on her qubits is 
one of 00, 01, 10 and 11, depending on her measure- 
ment outcome Bob's qubit will be one of the above 
four possible states. Once Alice sends to Bob her 
two-classical-bit measurement outcome through 
a classical channel, he applies appropriate quantum 
gates to his state and recovers \^p). For example, if 
her measurement is 00, Bob's state is 1-0), and he 
does not need to apply any quantum gate. If her 
measurement is 01, then Bob needs to apply a cr^ 
gate to his state ao|l) -|- ai|0) and yields {ip). If her 
measurement is 10, then applying a cr^ gate to his 
state ao|0) — ai|l) Bob recovers \^p). If her measure- 
ment is 11, then Bob can fix up his state ao|l) ~ 
ai|0) to recover \tp) by applying first a ctj. gate and 
then a cr^ gate. Here the and cr^ gates are de- 



10 



Y. WANG 



fined by Pauli matrices <t^. and <t^ given by (7). 
In summary, according to AHce's measurement out- 
come, applying some appropriate quantum gates to 
his qubit Bob will recover the state I'i/')- 

A few important remarks about quantum telepor- 
tation are in the line. First, quantum teleportation 
does not involve any transfer of matter or energy. 
Alice's particle has not been physically moved to 
Bob; only its state has been transferred. Second, af- 
ter the teleportation Bob's qubit will be on the tele- 
ported state, while Alice's qubit will become some 
undefined part of an entangled state. In other words, 
what the teleportation does is that a qubit was de- 
stroyed in one place but instantaneously resurrected 
in another. Teleportation does not copy any qubits, 
and hence is consistent with the no-cloning theo- 
rem (which forbids the creation of identical copies of 
an arbitrary unknown quantum state; see Wootters 
and Zurek (1982)). Third, in order to teleportate 
a qubit, Alice has to inform Bob of her measure- 
ment by sending him two classical bits of informa- 
tion. These two classical bits do not carry complete 
information about the qubit being teleported. If the 
two bits are intercepted by an eavesdropper, he or 
she may know exactly what Bob needs to do in or- 
der to recover the desired state. However, this infor- 
mation is useless if the eavesdropper cannot interact 
with the entangled particle in Bob's possession. Also 
the requirement of sending two bits of information 
via classical channel prevents quantum teleportation 
from transmitting information faster than the speed 
of light. 

3.3.3 Bell's inequality The Bell test experiments 
are designed to investigate the validity of the en- 
tanglement effect in quantum mechanics through 
Bell's inequality. Over the past four decades many 
physical experiments on quantum systems were con- 
ducted to check the validity of Bell's inequality and 
resulted in some violation of the inequality. For ex- 
ample. Aspect, Grangier and Roger (1981, 1982a, 
1982b) provided overwhelming support to the vi- 
olation of Bell's inequality. The experimental re- 
sults are often invoked as the proof of quantum non- 
locality and lack of realism that no particle has defi- 
nite form until it is measured and measuring a quan- 
tum entity can instantaneously influence another far 
away. See Aspect, Grangier and Roger (1981, 1982a, 
1982b), Bohm (1951), Bell (1964), Clauser et al. 
(1969) and Einstein, Podolsky and Rosen (1935). 
Below we describe the CHSH version of the Bell's 
inequality (Clauser et al. (1969)). 



Suppose Xi, z = 1,2,3,4, are four random vari- 
ables taking values ±1. Consider an ordinary exper- 
iment with two people, Alice and Bob. In the ex- 
periment Alice observes Xi or X2 while Bob meas- 
ures X3 or X4. Consider the quantity X1X3 + X2X3 + 
X2X4 — X1X4. It is equal to 

{Xi + X2)X3 + {X2 - Xi)X4 = ±2 < 2. 

Regardless of the distributions of Xi, taking expec- 
tation on both sides of the above inequality we arrive 
at the famous Bell inequality. 



(9) 



E{XiX3) + ^(^2X3) + E{X2X4) 
- E{XiXi) < 2. 



The violation of Bell's inequality demonstrates en- 
tanglement effect in quantum mechanics. In fact, 
quantum experiments yield a quantum version of 
the inequality. Consider that a quantum system of 
two qubits is prepared in a Bell state 

Alice takes the first qubit of while Bob gets its 
second qubit. Define four observables with eigenval- 
ues ±1, 

on the first qubit and 

— 7= , ^4 - 



V2 ' " V2 ' 

on the second qubit, where cTr^ and cr^ are Pauli ma- 
trices given by (7). Again Alice performs measure- 
ments on Xi or X2 while Bob measures X3 or X4. 
The quantum expectations of X1X3, X2X3, X2X4, 
X1X4 in the state I?/') are calculated below: 

E^(XiX3) = 4s, E^(X2X3) ^ 



V2' 



V2' 



E^(X2X4) = — ^, E^(XiX4) = — — ^. 

Here the observable product is in the sense of tensor 
product. Thus we obtain a value in the quantum 
framework for the analog quantity on the left-hand 
side of the Bell's inequality (9) 

E^(XiX3) +E^(X2X3) +Ev,(X2X4) 

-E^(XiX4) = 2V2, 

which exceeds 2 and hence violates the Bell's in- 
equality. In fact, the quantum version of the Bell's 
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inequality is the Tsirelson's inequality (Tsirelson, 
1980) which shows that in any quantum state p, 

Ep(XiX3) + Ep(X2X3) + Ep(X2X4) 

(10) 

- Ep(XiX4) < 2V2. 
3.4 Quantum Parallelism 

Quantum computation has an amazing feature 
termed as quantum parallelism, which may be heu- 
ristically explained by the following oversimplify- 
ing description: a quantum computer can simulta- 
neously evaluate the whole range of a function f{x) 
at many different values of x. 

For function f{x) with b bit input x = xi - ■ ■ Xj, 
and 1 bit output f{x), we illustrate quantum par- 
allel evaluation of its values at many different x si- 
multaneously as follows. First we apply b Hadamard 
gates to the first b qubits of |0 • • • 0)|0) to obtain 



|0) + |1) |0) + |1) 



V2 



V2 



|0) 



1 



|x)|0), 



X — X'^ ' ' ' Xf~f^ j — 

where the sum is over all possible 2^ values of x. 
Second, apply quantum circuit U/ defined in (6) to 
the obtained 6-1-1 qubit state to yield 



The quantum circuit with b Hadamard gates is ex- 
tremely efficient in producing an equal superposi- 
tion of all 2* computational basis states with only b 
gates; and quantum parallelism enables simultane- 
ous evaluation of the whole range of the function /, 
although we evidently evaluate / just once with sin- 
gle quantum circuit U f applied to the superposition 
state. To make it more clear we consider the case of 
6 = 1. Apply circuit U/ to a superposition state as 
follows: 

u/M^io>^-T,yi°°)+i'°)' 



|0/(0)) + |l/(l)) 
V2 



One application of a single circuit U j results in a su- 
perposition state whose two components contain in- 
formation about both /(O) and /(I), as if we have 
evaluated f{x) at values and 1 simultaneously. 
The quantum parallelism is in contrast with clas- 
sical parallelism, where multiple circuits each built 
to compute one value of f{x) are executed simulta- 
neously. Quantum parallelism arises from superpo- 
sition states. A superposition state has many com- 



ponents, each of which may be thought of as a sin- 
gle argument to function f{x). Because of quantum 
nature, a single circuit Uj applied once to the su- 
perposition state is actually performed on each of 
the components of the superposition, and the whole 
range of the values of function f{x) is stored in the 
resulted outcome superposition state. 

The quantum parallelism can be a potentially pow- 
erful tool for computational statistics. For example, 
Bayesian analysis often encounters the problems of 
evaluating sums over 2^ quantities, with 6 propor- 
tional to sample size or the number of variables. For 
moderate to large 6, the evaluation of such sums is 
computationally prohibitive by classical computers 
(Vidakovic (1999)). Because of the quantum paral- 
lelism, it is possible for quantum computers to per- 
form such computing tasks. 

4. QUANTUM ALGORITHMS 

Quantum algorithms are described by quantum 
circuits that take input qubits and yield output mea- 
surements for the solutions of the given problems. 
As a classical algorithm is a step-by-step problem- 
solving procedure, with each step performed on a clas- 
sical computer, a quantum algorithm is a step-by- 
step procedure to solve a problem, with each step ex- 
ecuted by a quantum computer. Although all classi- 
cal algorithms can also be carried out on a quantum 
computer, we refer to quantum algorithms as the 
algorithms that utilize essential quantum features 
such as quantum superposition and quantum entan- 
glement. While it is true that all problems solvable 
on a quantum computer are solvable on a classi- 
cal computer, and problems undecidable by classi- 
cal computers remain undecidable on quantum com- 
puters, what makes quantum algorithms exciting is 
the faster speed that quantum algorithms might be 
able to achieve, compared to classical algorithms, 
for solving some tough problems. The well-known 
quantum algorithms are Shor's factoring algorithm 
and Grover's search algorithm. Shor's algorithm and 
Grover's algorithm run, respectively, exponentially 
faster and quadratically faster than the best known 
classical algorithms for the same tasks. Common 
techniques used in quantum algorithms include quan- 
tum Fourier transform, phase estimation and quan- 
tum walk. 

4.1 Quantum Fourier Transform 

The quantum Fourier transform is defined to be 
a linear transformation on n qubits that maps the 
computational basis states | j), j = 0, 1, . . . , 2" — 1, 
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to superposition states as follows: 



2"-l 
on Z — / 



ijk/2" I 



A:=0 



The inverse of quantum Fourier transform is given 
by 



1 



/on Z— / 

^ i=0 



|J>- 



We use the binary representation to express the 
state j = jiT"-^ + j22'^-2 + • • • + j„20 as j = ji • • • i„ 
and represent binary fraction + j£_|_i/2^ + • • • + 
j^/2'"-^+i as O.j^j^+i • • • j™, where 1 < £ < m < 2". 
Then the quantum Fourier transform of state |j) = 
|jij2'"Jn) has the following useful product repre- 
sentation: 

|jii2---Jn)^^(|0)+e2-0->'|l)) 

•(|0)+e2'^*°-^"-i^"|l))--- 

.(|0)+e2'^*°-J'iJ2-^"|l)). 

It can be easily checked from the product represen- 
tation that with quantum parallelism the quantum 
Fourier transform can be realized as a quantum cir- 
cuit with only O(n^) operations, while classically the 
fast Fourier transform requires 0(n2") operations 
for processing 2"- data, which indicates an exponen- 
tial speed-up (Nielsen and Chuang (2000)). Realiz- 
ing such an exponential saving accommodated by 
quantum parallelism requires clever measurement 
schemes. Successful examples include quantum phase 
estimation and Shor's algorithms for factoring and 
discrete logarithm. 

4.2 Phase Estimation 

Quantum algorithms are of random nature in the 
sense that they are able to produce correct answers 
only with some probabilities. Consider quantum 
phase estimation which provides the key to many 
quantum algorithms. Assume that a unitary opera- 
tor U has an eigenvector \x) with eigenvalue e^"^^^ . 
The phase ip of the eigenvalue is unknown and the 
goal of the phase estimation algorithm is to esti- 
mate if based on the assumption that the state |x) 
can be prepared and the controlled-U^^ operations 
[see Section 3.2 for control gate] can be performed 
for suitable nonnegative integers j. 

The registers are used in phase estimation. The 
first register consists of h qubits initially in the state 
|0). The second register starts in the state |x) and 



involves enough qubits to store \x). The phase esti- 
mation procedure is performed in two stages. 

First, we apply Hadamard transform to the first 
register and then controlled-U operations on the 
second register, with U raised to successive pow- 
ers of 2, to obtain the final state with the second 
register unchanged and the first register given by 

1 



2^2 



(|0) + e- 



27142" 



(11) 



•(|0) + e 



27142" 



^ll)) 

'^\l))---{\Q) + e^^''^°^\l)) 



1 



fc=0 



If is expressed exactly in h bits as ip = G.^pi ■ ■ ■ ipb, 
(11) becomes 



1 



2V2 



(|0) +e2'^^0-'^<'|l)) 



(12) 



.(|0)+e2^^'°-^"-i^"|l))... 

• (|0) +e2^*°-^i^i-^"|l)), 

which is the quantum Fourier transform of the prod- 
uct state \fiip2---^b)- 

The second stage of phase estimation is to take 
the inverse quantum Fourier transform on the first 
register. For ip = O.ipi • • • ipi,, the output state from 
the second stage is \^piip2 - • ■ Pb)-, and a measurement 
in the computational basis yields ipi - ■ -ipi, and di- 
viding the measurement by 2^ gives cpi ■ ■ ■ = 
O.ifi ■ • ■ iph = if. We obtain a perfect estimate of tp. 

Now we consider the case that p cannot be ex- 
pressed exactly with a b bit binary expansion. Take 
< rj < 2^ to be the integer that its binary frac- 
tion r]/2^ = 0.?7ir/2 • • • % is the first b bit represen- 
tation in the binary expansion of p, which satisfies 
0<p-'i]/2^ <2-K 

Perform the inverse quantum Fourier transform 
on the first register given by (11), which is obtained 
in the first stage results, and get 



^ 2°-l 



2f' 



^-2iTik£/2'' ^2iTi<pk 



kl=0 



2"-l 

E 

1=0 



where amplitudes of [(ry -|- ^)(mod2^)) are 



2°-l 



2^ 
1 



k=0 



27ri[(/j-(r)+£)/2"]|A: 



_ ^2TTi{2''^-ri-e) 
o27ri(v3-rj/2"-£/2") 
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Assume that the result of the final measurement 
from phase estimation is fj and dividing the result by 
2* gives = 77/2*. Let ( be the specified accuracy for 
the phase estimation procedure. By adding up |/3^p 
with i being within (^2'', we bound the probability 
that the obtained (f is within from ip: 

P{\^-ip\<C)>Pi\v-v\<C2''-i) 

1 



> 1 



2(C2'' - 2) 











+ 


log, (2 + 1) 



For e > 0, set 



(13) 



Then P{\(p — ip\ < C) ^ 1 — that is, with probability 
at least 1 — e the phase estimation procedure can 
successfully produce (f within from the true ip. 
See Nielsen and Chuang (2000). 

4.3 Statistical Analysis 

The phase estimation algorithm requires b qubits 
for the first register to achieve [— log2 C] bit accuracy 
and success probability 1 — e. With accuracy fixed, 
to increase the success probability the required qu- 
bits 

6~ 1 - log2C+ . } 

4e log 2 

which grows at a very fast rate. For example, an 
increase in success probability from 90% to 99% re- 
quires eighteen times of qubit increase compared to 
the change from 80% to 90%. 

Quantum algorithms are of random nature in the 
sense that they often produce correct answers only 
with certain probabilities. The success probabilities 
depend upon the schemes of the algorithms as well 
as the context of applications. Given a quantum al- 
gorithm for solving a problem, a common practice is 
to repeatedly run the quantum algorithm to achieve 
high probability of successfully obtaining a correct 
answer. Consider that the phase estimation proce- 
dure is repeatedly run n times to obtain results 
(fii, . . . ,ipn- Then (fi, . . . ,(fn may be treated as i.i.d. 
random variables with each ipj satisfying 

We may statistically model tpj by the gross error 
model (Huber and Ronchetti (2009)) as follows. As- 
sume that tpj are independently and identically gen- 
erated from (1 — e)F{x) + eH{x), where F{x) is the 
distribution of the correct answers that are within (" 
from true and H{x) is the distribution of wrong 
answers that are at least (" away from true f. Then (pj 



are correct with probability 1 — e and incorrect with 
probability e. 

If the outcome result of the algorithm is verifi- 
able to be a correct answer or not [as in the case 
of Shor's algorithms for factoring and order-finding 
in Section 4.4 below], the obtained result from each 
run is checked to be a correct answer or not. Then 
the number of times required to run the algorithm 
for obtaining a correct answer follows a geometric 
distribution. Thus the probability that we obtain 
a correct answer in n repetitions is equal to 

P (obtain a correct answer in n trials) 

= 1 — i-*(no success in the n trials) = 1 — e". 

Since e" goes to zero geometrically fast, we may 
choose a moderate e with fewer qubits to achieve 
very high probability of successfully obtaining a cor- 
rect answer by repeatedly running the algorithm 
enough times. 

On the other hand, if the outcome result is not 
verifiable to be a correct answer or not [as in the 
case of phase estimation], careful analysis is needed 
to design ways for obtaining a correct answer with 
very high probability. As wrong answers are far away 
from true estimators like sample average of (^1, ... , 
(fn may not estimate well. We adopt a robust sta- 
tistical method to estimate (p by a-trimmed mean ip, 
which is defined as follows. Ordering ipi, . . . ,(pn and 
then removing [na] largest and [na] smallest ones, 
we take the average of the remaining (pj as a-trimmed 
mean, where a is chosen to be greater than e/2. 
One example is the sample median of ipi, . . . ,(pn. 
The probability that p is within from ip can be 
calculated from the binomial probability as follows: 

P(|(^ — V^l < C) ^ P(more than n(l — 2a) number 

of (pj are within ( from p) 



n 
k 

fc=[n(l-2a)]-l 

As n — )■ 00, ne approaches to infinity. The binomial 
probability can be approximated by resorting to a nor- 
mal approximation, yielding 

n 



fc=[n(l-2a)]-l 



[l-ere 



k ^n—k 



where <!>(•) is the standard normal distribution func- 
tion. Since 2a — e > 0, as n increases, P{\(p — ip\ < Q 
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approaches to 1 exponentially fast. Combining the 
two cases together, we arrive at the following theo- 
rem. 

Theorem 1. Suppose that the outcome ip of 
a quantum algorithm obeys the gross error model 
that with probability 1 — e it produces a correct an- 
swer and probability e it gives a wrong answer. Then 
by repeatedly running the quantum algorithm we will 
obtain a correct answer with probability approaching 
1 exponentially fast in the number of repetitive runs. 

For a quantum algorithm that produces a correct 
answer with probability 70% and a = 0.2, in order 
to obtain a correct answer with 0.999 probability we 
need to run the algorithm five times and 20 times, 
respectively, for the cases that the outcome results 
are verifiable and not verifiable. 

4.4 Factoring and Order-Finding Algoritlims 

The factoring problem is to find all prime fac- 
tors of a given positive composite number such that 
the product of these prime numbers is equal to the 
composite number. Factoring is known to be a very 
hard problem for classical computers. Shor (1994, 
1997) developed a quantum algorithm for the fac- 
toring problem that is exponentially faster than the 
most efficient known classical factoring algorithm. 

Shor's quantum algorithms work as follows. Math- 
ematically the factoring problem is equivalent to 
the order-finding problem that for two positive inte- 
gers X and N , X < N , with no common factors, find 
the smallest integer r such that dividing x^' by N 
we obtain a reminder 1 (Shor (1997); Nielsen and 
Chuang (2000)). The quantum algorithm for factor- 
ing is reduced to a quantum algorithm for order- 
finding. The quantum algorithm for order-finding is 
to apply the phase estimation algorithm to the uni- 
tary operator 

U|y) = \xy{modN)) . 
The eigenvectors of U are 

I \ 1 / —2Trisk\ i. 

I'"^/ = 7FZ^'^''pI r j'"^ modiV), 

^ k=o ^ ^ 

s = 0,1, . . . ,r — l,i = 

with corresponding eigenvalues exp(27rzs/r). Using 
the phase estimation algorithm we can obtain the 
eigenvalues exp(27rzs/r) with high accuracy and thus 
find the order r with certain probability. 

While the quantum factoring algorithm can ac- 
complish the task of factoring an n-bit integer with 



operations of order log n log log n, the current best 
known classical algorithm requires operations of or- 
der exp(n^/^ log^/^ n) to factor an n-bit composite 
number (Crandall and Pomerance (2001)). Note that 
the number of operations required in the best classi- 
cal algorithm grows exponentially in the size of the 
number being factored. Because of the exponential 
complexity, the factoring problem is generally re- 
garded as an intractable problem on classical com- 
puters. 

The factoring problem plays an important role in 
cryptography. Cryptography is to enable two par- 
ties, Alice and Bob, to communicate privately, while 
it is very difficult for the third parties to "eaves- 
drop" on the contents of the communications. Ex- 
amples include ATM cards, computer passwords, in- 
ternet commences, clandestine meetings and mili- 
tary communications. Two cryptographic protocols 
used in the communications are private key cryp- 
tosystem and public key cryptosystem. A private 
key cryptosystem requires the two communicating 
parties to share a private key. Alice uses the key 
to encrypt the information, sends the encrypted in- 
formation to Bob who uses the key to decrypt the 
received information. The severe drawback of the 
private key cryptosystem is that the parties have 
to safeguard the key transmission from being eaves- 
dropped. A public key cryptosystem invented in the 
1970s requires no sharing secret key in advance. Bob 
publishes a "public key" available to the general 
public, and Alice uses the public key to encrypt in- 
formation and sends the encrypted information to 
Bob. The encryption transformation is specially cre- 
ated such that with only the public key, it is ex- 
tremely difficult, though not impossible, to invert 
the encryption transformation. When publishing the 
public key Bob keeps a matched secret key for easy 
inversion of the encryption transformation and de- 
cryption of the received information. One of the 
most widely used public key cryptosystems is the 
RSA cryptosystem, which is named after its cre- 
ator Rivest, Shamir and Adleman (Menezes, van 
Oorschot and Vanstone (1996); Rivest, Shamir and 
Adleman (1978)). RSA is built on the mathematical 
asymmetry of factoring: it is easy to multiply large 
prime numbers and obtain their product as a com- 
posite number but hard to find the prime factors 
of a given large composite number. RSA encryption 
keeps the large primes as a secret key and uses their 
product to make a "public key." Because of its expo- 
nential complexity, tremendous efforts tried to break 
the RSA system so far have resulted in vain, and 
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there is a widespread belief that the RSA system is 
secure against any classical computer based attacks. 
As the factoring problem can be efficiently solved 
by Shor's quantum factoring algorithm, a quantum 
computer can break the RSA system easily. Fortu- 
nately, while quantum mechanics takes away with 
one hand, it gives back with the other. A quan- 
tum procedure known as quantum cryptography or 
quantum key distribution can do key distribution so 
that the communication security cannot be compro- 
mised. The idea is based on the quantum principle 
that observing a quantum system will disturb the 
system being observed. If there is an eavesdropper 
during the transmission of the quantum key between 
Alice and Bob, eavesdropping will disturb the quan- 
tum communication channel that is used to estab- 
lish the key, and the disturbance will make eaves- 
dropping visible. Alice and Bob will throw away the 
compromised key and keep only the secured key for 
their communication. 

4.5 Quantum Search Algorithm 

Suppose that you would like to find the name 
corresponding to a given phone number in a tele- 
phone directory; or suppose that there are some lo- 
cations in a given city you would like to visit and 
wish to find the shortest route passing through all 
the locations. If there are names in the telephone 
directory or possible routes to pass through all 
the locations, search algorithms by classical com- 
puters usually require operations of order A^. One 
such simple classical algorithm is to check exhaus- 
tively all names to find a name matching with the 
given phone number or to search all possible routes 
and then find the shortest route among all routes. 
However, Grover (1996, 1997) developed a quantum 
search algorithm that needs only operations of order 
^/N to find a solution to the search problem. 

The quantum search algorithm works as follows. 
Suppose that the search space has A^ elements and 
the search problem has exactly M solutions. Assume 
M < N/2. (For the silly case of M > N/2, we ei- 
ther search for the solution by doing random selec- 
tion from the search space or double the number of 
the elements in the search space by adding A^ extra 
non-solution elements to the search space.) The al- 
gorithm works by creating superposition state with 
Hadamard gate, 

x=0 



and then applying a so-called Grover iteration (or 
operator) repeatedly. Set 

where the summations over x' and x" denote sums 
over all non-solutions and solutions, respectively. 
Then we can express \'^) as follows: 

, IN-M.^. [m. . 

The Grover operator is to perform two reflections, 
one about the vector and another about the vec- 
tor \tp). The two reflections together are a rotation 
with angle 9 in the two-dimensional space spanned 
by \4>) and |<^), where 

In -M 

cos(0/2) = y^^. 

After the rotation, the initial state = cos (0/2)10) -|- 
sm{6 /2)\ip) becomes state 

cos(36'/2)|(^) -Fsin(36l/2)|v9). 

Thus each application of the Grover operator is a ro- 
tation with angle 9. The initial state IV') has angle 
7r/2 — 9/2 with \{p)] after the first rotation, the re- 
sulted state has angle 7r/2 — 30/2 with \(p)] and in 
general after the rth rotation, the resulted state has 
angle 7r/2 — (2r -|- 1)9/2 with \ip). Repeatedly apply- 
ing the Grover operator, we rotate the state vector 
near \(p). With the initial state = cos{9 /2)\(f)) + 
s\n{9 /2)\ip) , we need to rotate through arccos \fMjN 
radians to transform the state vector to After 
R = arccos (v^M/A^)/^ = 0{^N/M) times of appli- 
cations of the Grover operator, we rotate the state 
vector 1-0) to within an angle 9/2 of \(p). Performing 
measurements of the state yields a solution to the 
search problem with probability at least cos^(0/2) > 
1 - M/N. 

The number of iterations R depends on M, the 
number of solutions. Since R < ■k/{20) and 0/2 > 
sin(0/2) = ^/^TJN, R < {n/A)./N/M. Typically, 
M < A^, f« sin « 2^M/N, thus i? (vr /4) • ^N/M. 
We estimate the number of solutions by quantum 
counting, which is to combine the Grover operator 
with the phase estimation method. Under the basis 
10) and |(/?) the Grover operator has eigenvalues e*^ 
and e*(^'^~^^ . Applying the phase estimation method 
we can estimate the eigenvalues and thus 9 with 
prescribed precision and probability, which in turn 
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yields M. The combination of the quantum count- 
ing and search procedure will find a solution of the 
search problem with certain probability. Repeating 
the quantum search algorithm will boost the prob- 
ability and enable us to obtain a solution to the 
search problem. 

Quantum walk and quantum Markov chain are 
currently being investigated for new quantum search 
algorithms and quantum speed-up of Markov chain 
based algorithms (Aharonov and Ta-Shma (2003), 
Childs et al. (2003); Childs (2010); Tulsi (2008); 
Shenvi, Kempe and Whaley (2003) and Szegedy, 
2004). In Section 5 we show that the quantum search 
algorithm can also be viewed as a quantum simula- 
tion procedure. 

5. QUANTUM SIMULATION 

Quantum simulation is to intentionally and artifi- 
cially mimic a natural quantum dynamics, which is 
hard to access, and analyze, by a computer-gener- 
ated quantum system, which is easy to manipulate 
and investigate. It provides scientific means for sim- 
ulating complex biological, chemical or physical sys- 
tems in order to study and understand certain scien- 
tific phenomena and evaluate hard-to-obtain quan- 
tities in the systems. Examples in modern scien- 
tific studies include the estimation of dielectric con- 
stant, proton mass, and precise energy of molecu- 
lar hydrogen, the study of superconductivity, the 
test of novel nano- materials, and the design of new 
biomolecules. 

To simulate a quantum system we need to solve 
the Schrodinger equation (1) which governs the dy- 
namic evolution of the system. For a typical Hamil- 
tonian with real particles the Schrodinger equation 
usually consists of elliptical differential equations, 
each of which can be easily simulated by a classical 
computer. However, the real challenge in simulating 
a quantum system is to solve the exponential num- 
ber of such differential equations. For a quantum 
system of h qubits, its states have 2^ amplitudes. To 
simulate the dynamic behavior of h qubits evolving 
according to the Schrodinger equation, we need to 
solve a system of 2^ differential equations. Due to 
the exponential growth in the number of differential 
equations, the simulation of general quantum sys- 
tems by classical computers is very inefficient. Clas- 
sical simulation of quantum systems is feasible for 
the cases where insightful approximations are avail- 
able to dramatically reduce the effective number of 
equations involved. Quantum computers may excel 
in simulating physically interesting and important 



quantum systems for which efficient simulation by 
classical computers may not be available. 

5.1 Simulate a Quantum System 

The key of quantum simulation is to solve the 
Schrodinger equation (1) which has solution 

(14) m))=e-'^'\m), i=^- 

Numerical evaluation of e~*^* is needed. The Hamil- 
tonian H is usually exponentially large and extremely 
difficult to exponentiate. The common approach in 
numerical analysis is to use the first-order linear ap- 
proximation, 1 — iH(5, of e"*^*^*"'"'') — e~*^*, which 
often yields unsatisfactory numerical solutions. 

Many classes of Hamiltonians have sparse repre- 
sentations. For such sparse Hamiltonians we can find 
efficient evaluation of the solutions (14) with high- 
order approximation. For example, the Hamiltoni- 
ans in many physical systems involve only local in- 
teractions, which originate from the fact that most 
interactions fall off with increasing distance in lo- 
cation or increasing difference in energy. In the lo- 
cal Hamiltonian case, the Hamiltonian of a quantum 
system with a particles in a ci-dimensional space has 
the form 

L 

(15) H = 2j^H,, 

1=1 

where L is a polynomial in a + d, and each acts 
on a small subsystem of size free from a and d. For 
example, the terms are typically two-body inter- 
actions and one-body Hamiltonians. Hence e~^^^^ 
are easy to approximate, although e~*^^ is very hard 
to evaluate. Since and are non-commuting, 
g-jH(5 _^ g-iHi5 . . . g-jHL<5^ Applying a modification 

of the Trotter formula (Kato (1978); Trotter (1959); 
Yu (2001)) we obtain 

e-''^^ = Us + 0{5^), 

where 

(16) Us = [e^^^i-^ • • • e"^"^'^] [e"*"^"^ • • • e"''^'^] . 

Thus we can approximate e~^^^ by which needs 
to evaluate only each e~^^^^ . 

Assume that the quantum system starts at t = 
with initial state |V'(0)) and ends at final time t = 1. 
For an integer m, set tj = j/m, j = 0, 1, . . . , m. The 
quantum simulation is to apply approximation of 
g-2iH<5 evaluate (14) at tj iteratively and generate 
approximate solutions to \ip{tj)). Denote by \ip{tj)) 
the state at tj obtained from the quantum simula- 
tion as an approximation of the true state \'ip{tj)) 
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at tj. Then for j = 1, . . . , m, 



(17) 



-2iH<5| 



(i.-i)) 



mto)), 



While classical computers are inefficient in simu- 
lating general quantum systems, quantum comput- 
ers can efficiently carry out the quantum simulation 
procedure and provide an exponential speedup for 
the quantum simulation on classical computers. In 
spite of the inefficiency, classical computers are cur- 
rently being used to simulate quantum systems in 
biochemistry and material science. Quantum simu- 
lation will be among the important applications of 
quantum computers. See Abrams and Lloyd (1997), 
Aspuru-Guzik, Dutoi and Head-Gordon (2005), Ben- 
nett et al. (2002), Berry et al. (2007), Freedman, Ki- 
taev and Wang (2002), Jane et al. (2003), Boghosian 
and Taylor (1998), Lloyd (1996), Nielsen and Chuang 
(2000), and Zalka (1998). 

5.2 Recast Quantum Search Algorithm as 
Quantum Simulation 

Grover's search algorithm discussed in Section 4.5 
is an important finding in quantum computation. It 
can be heuristically sketched quantum simu- 
lation by writing down an explicit Hamiltonian H 
such that a quantum system evolves from its initial 
state \tp) to \x) after some specified time, where x 
is a solution of the search problem. Of course the 
Hamiltonian H depends on the initial state {ip) and 
solution X. Suppose that \y) is another state such 
that |x) and \y) form an computational basis, and 
\ip) = a\x) + for real a and /3 with + = 1. 
Define Hamiltonian 



H = |2;)(x| + = I + a{l3a^ + aa,), 

where cTx and cr^ are Pauli matrices defined in (7). 
Then 

exp(-im)lV') 

= e~**[cos(at)|V') — ism{at){Pax + acr^)!^)] 

= e~**[cos(at)|^) — ism{at)\x)]. 

Measuring the system at time t = 7r/(2a) yields the 
solution state \x). 

5.3 Quantum Monte Carlo Simulation 

Quantum theory is intrinsically stochastic and 
quantum measurement outcome is random. As many 
naturally occurring quantum systems involve a large 
number of interacting particles, due to the compu- 
tational complexity we are forced to utilize Monte 



Carlo techniques in the simulations of such quantum 
systems. The combination of Monte Carlo methods 
with quantum simulation makes it possible to ob- 
tain reliable quantifications of quantum phenomena 
and estimates of quantum quantities. Such combi- 
nation procedures are often referred to as quantum 
Monte Carlo simulation (Nightingale and Umrigar 
(1999); Rousseau (2008)). Consider the problem of 
estimating the following quantity: 



(18) 



e = Tr{Xp) = EiX), 



where X is an observable, X is its measurement re- 
sult, and p is the state of the quantum system under 
which we perform the measurements and evaluate 
the quantity 9. 

As p is the true final state of the quantum system, 
we denote by p the final state of the quantum sys- 
tem obtained via quantum simulation. The quantum 
systems are prepared in initial state Pq, and we use 
the quantum simulation procedure described above 
to simulate the evolutions of the systems from initial 
state pQ to final state p according to Schrodinger's 
equation (14) with Hamiltonian H given by (15). 
We repeatedly perform the measurements of such n 
identically simulated quantum systems at the final 
state and obtain measurement results Xi, 
We estimate 9 defined in (18) by 



(19) 



1 " 

n ^ ■' 



The target 9 given by (18) is defined under the 
true state p, while the simulated quantum system is 
under approximate final state p which is close to p. 
The measurement results Xi , . . . , X^ are obtained 
via quantum simulation from the quantum systems 
in the simulated state p. Therefore, the Monte Carlo 
quantum estimator 9 in (19) involves both bias and 
variance. Wang (2011) studied the quantum simu- 
lation procedure and investigated the bias and vari- 
ance of 9. The derived bias and variance results can 
be used to design optimal strategy for the best uti- 
lization of computational resources to obtain the 
quantum Monte Carlo estimator. 

6. QUANTUM INFORMATION 

Classical information theory is centered on Shan- 
non's two coding theorems on noiseless and noisy 
channels. The noiseless channel coding theorem quan- 
tifies the number of classical bits required to store 
information for transmission by Shannon entropy, 
while the noisy channel coding theorem quantifies 
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the amount of information that can be rehably trans- 
mitted through a noisy channel by an error-correction 
coding scheme. The quantum analogs of Shannon 
entropy and Shannon noiseless coding theorem are 
von Neumann entropy and Schumacher's noiseless 
channel coding theorem, respectively. The von Neu- 
mann entropy is defined to be S{p) = — tr(/9logp). 
Schumacher's noiseless channel coding theorem quan- 
tifies quantum resources required to compress quan- 
tum states by von Neumann entropy (Schumacher 
(1995)). Analogous to Shannon's noisy channel cod- 
ing theorem, a theorem known as Holevo-Schuma- 
cher-Westmoreland theorem can be used to com- 
pute the product quantum state capacity for some 
noisy channels (Holevo (1998); Schumacher and West- 
moreland (1997)). However, communications over 
noisy quantum channels are much less understood 
than the classical counterpart. It is an unsolved prob- 
lem to determine quantum channel capacity or the 
amount of quantum information that can be reli- 
ably transmitted over noisy quantum channels. See 
Hayashi (2006) and Nielsen and Chuang (2000). 

In spite of the above similarity, there are intrinsic 
differences between classical information and quan- 
tum information. Classical information can be dis- 
tinguished and copied. For example, we can identify 
different letters and produce an identical version of 
a digital image for back-up. However, quantum me- 
chanics does not allow unknown quantum states to 
be distinguished or copied exactly. For example, we 
cannot reliably distinguish between quantum states 
|0) and (|0) -|- |l))/\/2. If we perform measurement 
for quantum state |0), the measurement result will 
be with probability 1, while measuring quantum 
state (|0) -|- |l))/\/2 yields measurements or 1 with 
equal probability. A measurement result of can- 
not tell the identity of the quantum state being 
measured. A theorem known as a no-cloning theo- 
rem states that unknown quantum states cannot be 
copied exactly (Wootters and Zurek (1982); Nielsen 
and Chuang (2000)). 

As we discussed in Section 3.3, quantum entan- 
glement plays a crucial role in strange quantum ef- 
fects such as quantum teleportation, violation of 
Bell's inequality, and superdense coding (Hayashi 
(2006); Nielsen and Chuang (2000)). Entanglement 
is a new type of resource that differs vastly from 
the traditional resources in classical information the- 
ory. We are far from having a general theory to 
understand quantum entanglement but encouraging 
progress made so far reveals the amazing property 



and intriguing structure of entangled states and re- 
markable connections between noisy quantum chan- 
nels and entanglement transformation. Consider 
quantum error-correction for reliable quantum com- 
putation and quantum information processing. Quan- 
tum error-correction is employed in quantum com- 
putation and quantum communication to protect 
quantum information from loss due to quantum noise 
and other errors like faulty quantum gates. Classi- 
cal information uses redundancy to achieve error- 
correction, but the no-cloning theorem presents an 
obstacle to copying quantum information and for- 
mulating a theory of quantum error-correction based 
on simple redundancy. Again quantum entanglement 
comes to the rescue. It is forbidden to copy qubits 
but we can spread the information of one qubit onto 
a highly entangled state of several qubits. Shor (1995) 
first discovered the method of formulating a quan- 
tum error-correction code by storing the information 
of one qubit onto a highly entangled state of nine 
qubits. Over time several quantum error-correction 
codes are proposed (Calderbank and Shor (1996); 
Cory et al. (1998); Steane, 1996a, 1996b). These 
quantum error-correction codes can protect quan- 
tum information against quantum noise, and thus 
quantum noise likely poses no fundamental barrier 
to the performance of large-scale quantum comput- 
ing and quantum information processing. 

Here is how quantum error-correction codes work. 
We consider the single qubit case. First assume that 
a qubit ao|0) +ai|l) is passed through a bit fiip 
channel which fiips the state of a qubit from |0) 
to |1) and from |1) to |0), each with probability p, 
and leaves each of states |0) and |1) untouched with 
probability 1 — p. We describe a bit flip code that 
protects the qubit against quantum noise from the 
bit flip channel. 

We encode states |0) and |1) in three qubits, with 
|0) encoded as |000) and |1) as |111). Thus the qubit 
state ao|0) -|- ai|l) is encoded in three qubits as 
ttolOOO) -I- ai|lll). We pass each of the three qubits 
through an independent copy of the bit fiip chan- 
nel, and assume that at most one qubit is flipped. 
The following simple two-step error-correction pro- 
cedure can be used to recover the correct quantum 
state. 

Step 1. Perform a measurement on a specially con- 
structed observable and call the measurement result 
an error syndrome. The error syndrome can inform 
us what error, if any, occurred on the quantum state. 
The observable has eigenvalues 0, 1, 2 and 3, with 
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corresponding projection operators, 

Qo = |000)(000| + |111)(111| no error, 
Qi = |100)(100| + |011)(011| 

bit flip on the first qubit, 
Q2 = |010)(010| + |101)(101| 

bit flip on the second qubit, 
Q3 = |001)(001| + |110)(110| 

bit flip on the third qubit. 

If one of three qubits has one or no bit flip, the error 
syndrome will be one of 0, 1, 2 and 3, with corre- 
sponding to no flip, and 1, 2 and 3 to a bit flip on 
the first, second and third qubit, respectively. For 
example, if the first qubit is flipped, the corrupted 
state is IV') = ao|100) + ai|011). Since (V'lQilV') = 1 
and (V'lQjIV') = for j 7^ 1, in this case the error 
syndrome is 1. Although performing measurements 
usually causes change to the quantum state, the 
speciality of the constructed observable is that syn- 
drome measurement does not perturb the quantum 
state: it is easy to check that the state is {ip) both 
before and after the syndrome measurement. While 
the syndrome provides information about what flip 
error has occurred, it does not contain any infor- 
mation about the state being protected, that is, it 
does not allow us to deduce anything about the am- 
plitudes ao and ai. Such a special property is the 
generic feature of syndrome measurement. 

Step 2. The error type supplied by the error syn- 
drome can inform us what procedure to use to re- 
cover the original state. For example, error syndrome 
1 indicates a bit flip on the first qubit, and a flip on 
the first qubit again will perfectly recover the orig- 
inal state ctolOOO) +ai|lll). The error syndrome 
implies no error and doing nothing, and error syn- 
dromes 1, 2 and 3 correspond to a bit flip again 
on the first, second, third qubit, respectively. The 
procedure will recover the original state with per- 
fect accuracy, if there is at most one bit flip in the 
encoded three qubits. The probability that more 
than one bit flipped is + 3p^(l — p) = 3p^ — 2p^, 
which is much smaller than the error probability p 
of making no-correction for the typical bit flip chan- 
nel. Thus the encoding and decoding scheme makes 
the storage and transmission of the qubit more reli- 
able. 

Next we consider a more interesting noisy quan- 
tum channel: a phase flip channel which, with proba- 
bility p, changes a qubit state ao|0) + ai|l) to ao|0) — 



aijl), and with probability 1 — p, leaves alone the 
qubit. The following scheme is to turn the phase 
flip channel into a bit flip channel. Let |+) = (|0) + 
|l))/\/2 and |-) = (|0) - |l))/\/2 be a qubit basis. 
The phase flip channel leaves alone states |+) and 
|— ) with probability 1— p and changes |+) to |— ) 
and vice versa with probability p. In other words, the 
phase flip channel with respect to the basis | +) and 
|— ) acts just like a bit flip channel with respect to 
the basis |0) and Thus we encode |0) as |+ + +) 

and |1) as | ) for protection against phase flip 

errors. The operations for encoding, error-detection 
and recovery are the same as for the bit flip channel 
but with respect to the |+) and |— ) basis instead of 
the |0) and |1) basis. 

Last we describe Shor error-correction code. It is 
a combination of the three-qubit phase flip and bit 
flip codes. First use the phase flip code to encode 
states |0) and |1) in three qubits, with |0) encoded as 

+ and |1) as | ); next, use the three-qubit 

bit flip code to encode each of these qubits, with |+) 
encoded as (|000) + |lll))/\/2 and | — ) encoded as 
(|000) - |111))/V2. The resulted nine-qubit code has 
codeworks as follows: 

|000) + |000) + |000) + 
^'^ V2 V2 ' 

|000) - |000) - |000) - 

' 72 72 72 ■ 

With the mixture of both phase flip and bit flip 
codes, the Shor error-correction code can protect 
against bit flip errors, phase flip errors, as well as 
a combined bit and phase flip errors on any single 
qubit. In fact it has been shown that this simple 
quantum error-correction code can protect against 
the effects of any completely arbitrary errors on a sin- 
gle qubit (Shor (1995)). 

7. CONCLUDING REMARKS 

Quantum information science gains enormous at- 
tention in computer science, mathematics, physi- 
cal sciences and engineering, and several interdisci- 
plinary subfields are developing under the umbrella 
of quantum information. This paper reviews quan- 
tum computation and quantum information from 
a statistical perspective. We introduce concepts like 
qubits, quantum gates and quantum circuits in quan- 
tum computation and discuss quantum entangle- 
ment, quantum parallelism and quantum error-cor- 
rection in quantum computation and quantum infor- 
mation. We present major quantum algorithms and 
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show their advantages over the available classical 
algorithms. We illustrate quantum simulation pro- 
cedure and Monte Carlo methods in quantum simu- 
lation. As classical computation and simulation are 
ubiquitous nowadays in statistics, we expect quan- 
tum computation and quantum simulation will have 
a paramount role to play in modern statistics. This 
paper exposes the topics to statisticians and encour- 
ages more statisticians to work in the fields. There 
are many statistical issues in theoretical research as 
well as experimental work in quantum computation, 
quantum simulation and quantum information. For 
example, as measurement data collected in quan- 
tum experiments require more and more sophisti- 
cated statistical methods for better estimation, sim- 
ulation and understanding, it is imperative to de- 
velop good quantum statistics methods and quan- 
tum simulation procedures and study interrelation- 
ship and mutual impact between quantum estima- 
tion and quantum simulation. Since quantum com- 
putation is intrinsically random, and quantum simu- 
lation employs Monte Carlo techniques, as we point 
out in Section 4.3 and Wang (2011), it is impor- 
tant to provide sound statistical methods for ana- 
lyzing quantum algorithms and quantum simulation 
in general and study high-order approximations to 
exponentiate Hamiltonians and the efficiency of the 
resulted quantum simulation procedures in particu- 
lar. On the other hand, quantum computation and 
quantum simulation have great potential to revo- 
lutionize computational statistics. Below are a few 
cases in point. 

1. The "random numbers" generated by classical 
computers are pseudo-random numbers in the 
sense that they are produced by deterministic 
procedures and can be exactly repeated and per- 
fectly predicted given the deterministic schemes 
and the initial seeds. On the contrary, superposi- 
tion states enable quantum computers to produce 
genuine random numbers. For example, measur- 
ing (|0) + 1 1))/\/2 yields and 1 with equal prob- 
ability. In general we generate 6-bit binary ran- 
dom numbers = 0,1 as follows. 
Apply b Hadamard gates to b qubits of 1 • • • 0) 
to obtain 

|0) + |1) |0) + |1) _ 1 >r^| , 

X X'[^ ' ' ' X^if j ^5 

where the sum is over all possible 2^ values of x, 
and then measure the obtained qubits and yield 



fe-bit binary random numbers x = xi - ■ ■ xi, with 
equal probability. Quantum theory guarantees that 
such random numbers are genuinely random. Thus 
quantum computers are able to generate genuine 
random numbers and perform true Monte Carlo 
simulation. It is exciting to design general quan- 
tum random number generator and study quan- 
tum Monte Carlo simulation. Perhaps we may 
need to re-examine Monte Carlo simulation stud- 
ies conducted by classical computers. 

2. It is interesting to investigate the potential of 
quantum computation and quantum simulation 
for computational statistics. We expect that quan- 
tum computers may be much faster than classical 
computers for computing some statistical prob- 
lems. Moreover, quantum computers may be able 
to carry out some computational statistical tasks 
that are prohibitive by classical computers. Spe- 
cific examples are as follows: (a) We may use the 
basic ideas of Grover's search algorithm to de- 
velop fast quantum algorithms for implementing 
some statistical procedures. For example, find- 
ing the median of a huge data set is to search for 
a numerical value that separates the top and bot- 
tom halves of the data, and quantum algorithms 
can offer quadratical speedup for calculating me- 
dian and trimmed mean, (b) With genuine ran- 
dom number generator and faster mean evalu- 
ation, quantum computers may offer significant 
advantages over classical computers for Monte 
Carlo integration. For example, Monte Carlo in- 
tegration in high dimensions may be exponen- 
tially or quadratically faster on quantum com- 
puters than on classical computers, (c) It might 
be possible for quantum computers to carry out 
some prohibiting statistical computing tasks like 
the Bayesian computation discussed in Section 3. 
Some preliminary research along these lines may 
be found in Nayak and Wu (1999) and Heinrich 
(2003). 

3. As quantum computation and quantum simula- 
tion are ideal for simulating interacting parti- 
cle systems like the Ising model, it is fascinating 
to explore the interplay between quantum simu- 
lation and Markov chain Monte Carlo method- 
ology and the quantum potential to speed up 
Markov chain based algorithms. In fact, it has 
been shown that quantum walk based algorithms 
can offer quadratical speedup for certain Markov 
chain based algorithms (Magniez et al. (2011); 
Richter (2007); Szegedy (2004), and Wocjan and 
Abeyesinghe (2008)). 
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Finally we point out that quantum computers are 
wonderful but it is difficult to build quantum com- 
puters with present technology. To build a quantum 
computer the physical apparatus must satisfy re- 
quirements that the quantum system realized qubits 
needs to be well isolated in order to retain its quan- 
tum properties and at the same time the quantum 
system has to be accessible so that the qubits can be 
manipulated to perform computations and measure 
output results. The two opposing requirements are 
determined by the strength of coupling of the quan- 
tum system to the external entities. The coupling 
causes quantum decoherence. Decoherence refers to 
the loss of coherence between the components of 
a quantum system or quantum superposition from 
the interaction of the quantum system with its en- 
vironment. It is very crucial but challenging to con- 
trol a quantum system of qubits and correct the ef- 
fects of decoherence in quantum computation and 
quantum information. Quantum computing has wit- 
nessed great advances in recent years, and quantum 
computers of a handful of qubits and basic quan- 
tum communication devices have been built in re- 
search laboratories (see Barz et al. (2012); Clarke 
and Wilhelm (2008); DiCarlo et al. (2009); John- 
son et al. (2011); Lee et al. (2011); Neumann et 
al. (2008)), but there are technological hurdles in 
the development of a quantum computer of large 
capacity. History shows that scientific innovations 
and technological surprises are a never-ending saga. 
It is anticipated that quantum computers with a 
few dozen of qubits will be built in near future. As 
we have discussed in Section 3.1, such a quantum 
computer has capacity of a classical supercomputer. 
We are very optimistic that someday quantum com- 
puters will be available for statisticians to crunch 
numbers. For the time being, instead of waiting in 
the sidelines for that to happen, statisticians should 
get into the field of play. It is time for us to dive 
into this frontier research and work with scientists 
and engineers to speed up the arrival of practical 
quantum computers. As a last note, in 2011 a Cana- 
dian company called D-Wave has sold the claimed 
first commercial quantum computer of 128 qubits 
to the Lockheed-Martin corporation, despite the D- 
Wave's quantum system being criticized as a black 
box. Large scale quantum computers may be years 
away, but quantum computing is already here as a 
scientific endeavor to provoke deep thoughts and in- 
tegrate profound questions in physics and computer 
science. 
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